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ABSTRACT 


We perform 2.5D hybrid simulations with massless fluid electrons and kinetic 
particle-in-cell ions to study the temporal evolution of ion temperatures, temper¬ 
ature anisotropies and velocity distribution functions in relation to the dissipation 
and turbulent evolution of a broad-band spectrum of parallel and obliquely propa¬ 
gating Alfven-cyclotron waves. The purpose of this paper is to study the relative 
role of parallel versus oblique Alfven-cyclotron waves in the observed heating 
and acceleration of minor ions in the fast solar wind. We consider collisionless 
homogeneous multi-species plasma, consisting of isothermal electrons, isotropic 
protons and a minor component of drifting a particles in a finite-/? fast stream 
near the Earth. The kinetic ions are modeled by initially isotropic Maxwellian 
velocity distribution functions, which develop non-thermal features and temper¬ 
ature anisotropies when a broad-band spectrum of low-frequency non-resonant, 
u < 0.34f2 p , Alfven-cyclotron waves is imposed at the beginning of the simula¬ 
tions. The initial plasma parameter values, such as ion density, temperatures 
and relative drift speeds, are supplied by fast solar wind observations made by 
the Wind spacecraft at 1AU. The imposed broad-band wave spectra is left-hand 
polarized and resembles Wind measurements of Alfvenic turbulence in the solar 
wind. The imposed magnetic field fluctuations for all cases are within the iner¬ 
tial range of the solar wind turbulence and have a Kraichnan-type spectral slope 
a = —3/2. We vary the propagation angle from 9 = 0° to 6 = 30° and 6 = 60°, 
and find that the minor ion heating is most efficient for the highly-oblique waves 
propagating at 60°, whereas the protons exhibit perpendicular cooling at all 
propagation angles. 

Subject headings: micro-turbulence — solar wind — wave-particle interactions — 



oblique AIC waves 
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Introduction 


The dissipation of fluid-scale fluctuations in collisionless plasmas can occur via a 
turbulent cascade, followed by different kinds of wave-particle interactions. The partitioning 
of energy between minor ions, protons and electrons and the efficiency of the heating 
depends on the characteristics of the waves, the wave vector direction and the anisotropy 
of the fiuctuations carrying energy at small scales. Anisotropic fluctuations here refer to 
different wave power in parallel and perpendicular direction with respect to the background 
magnetic field. Spacecraft observations of magnetic field fluctuations in the solar wind often 
provide ambiguous single-point measurements, which makes it hard to uniquely determine 
the type of the observed fluctuations or their angle of propagation. While parallel Alfven 
waves and Alfvenic turbulence are ubiquitous in the solar wind and their existence has 


been unambiguously proven throughout the heliosphere (Bruno & Carbone 2005), much 


less is known about obliquely propagating waves in the interplanetary medium. The most 
common electromagnetic fluctuations, which can interact with the particles at the ion scales 
are parallel and oblique ion-cyclotron (ICWs), as well as kinetic Alfven waves (KAWs). 
Such plasma waves at kinetic scales are very hard to detect and identify with a single 
spacecraft as their interpretation relies on preliminary assumptions. Yet, recent analysis 
of solar wind data from Helios, MESSENGER, STEREO, Wind and Cluster spacecraft 
show a compelling evidence for the existence of parallel and oblique ion-cyclotron waves 


ICWs from 0.3AU to 1 AU (Jian et al. 2009, 2010 2014), as well as KAWs (Sahraoui et al. 


2010 Podesta & Gary 2011| He et al. 2012a Salem et al. 2012) near the Earth. Jian et al. 


(2009, 2010) have shown that the ICWs are rather dispersive with a wide range of wave 


numbers and frequencies, and statistically significant fraction of events propagate at small 
angles with respect to the direction of the ambient magnetic field. Lately, reconstruction 
of the magnetic helicity based on Ulysses and STEREO/IMPACT/MAG data from the 
undisturbed solar wind suggests that the solar wind turbulence at the ion scales can be 
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reconstructed by a superposition of parallel and oblique ICWs/KAWs (Podesta & Gary 


2011 He et al. 2012b). What role each of the different types of waves plays in relation to 


wave-particle interactions, anisotropic heating and differential acceleration is still an open 
question, which remains to be solved. The purpose of this paper is to compare the ion 
heating by parallel and oblique Alfven-cyclotron waves and determine which propagation 
angles are most efficient in preferentially heating the minor ions within the considered 
low-frequency turbulent wave spectra. Our secondary goal is to examine the evolutionary 
path of particles in the configuration space set by the ion temperature anisotropies and their 


plasma fa see (Bale et al. 2009 Maruca et al. 2012 Bourouaine et al. 2011, 2013 Wicks 


et al. 2013). To address these problems, we perform 2.5D hybrid simulations with fluid 


electrons, kinetic protons and a minor component of differentially streaming a particles. We 
investigate the turbulent evolution, wave-scattering and dissipation of initial broad-band 
spectra of parallel and oblique Alfven-cyclotron waves in collisionless finite plasma f3 fast 
solar wind conditions. The simulation setup is motivated by the magnetic field power 
spectra and the plasma properties as in situ observed by Wind spacecraft near the Earth. 
We start within the stable middle point where Tj_,/Tj|j = 1 and the plasma fa = 0.33 
for both ion species, see definitions in the next section. This point is far from all linear 


plasma instabilities thresholds for both protons and the a particles (Maruca et al. 2012 


Bourouaine et al. 2013) and yet the protons are cooled while the minor ions are heated 


in perpendicular direction and the parallel temperature of both species increases. Such 
cooling/heating can not be described by the Vlasov linear instability theory predictions. It 
is a nonlinear process, resulting from the turbulent wave-wave couplings and consequent 
wave-particle interactions with the minor ions, whose proper treatment requires direct 
numerical simulations. 


The aim of this paper is to provide the first steps towards comparison and understanding 
of the efficiency of parallel versus oblique small-scale Alfvenic fluctuations for the preferential 
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heating of minor ions in the fast solar wind. The initial spectrum resembles Alfvenic 
turbulence within the inertial range with a Kraichnan-type power spectral slope a = —3/2. 
The spectral index for the fluctuations relates to in situ observations and has been calculated 
based on the Magnetic Field Instrument (MFI) magnetometer data on board of Wind 
spacecraft in the undisturbed solar wind at 1AU. The electron temperature and the ion 
properties are obtained respectively from the the Three-Dimensional Plasma and Energetic 
Particle Investigation (3DP) instrument and the Solar Wind Experiment (SWE) particle 
detector on board of Wind. The next section presents the particle data and the magnetic 
field measurements used in our model to construct a consistent initial state for the hybrid 
simulations. The description of the 2.5D hybrid model and the theoretical reconstruction 
of the initial wave spectrum are presented in Section 3 and Section 4. The results of the 
numerical simulations are revealed in Section 5. Discussion and concluding remarks are 
presented in Section 6. 


2. Initial state for the simulations — Wind data 


In order to relate our model to real solar wind conditions we have constrained the 
initial state for the simulations to resemble in situ measurements in a finite-/? fast stream 
near the Earth. The initial plasma parameters for the simulations are based on analysis of 
Wind spacecraft observations in the undisturbed fast solar wind at 1AU. In what follows 
below we will use the 92 seconds average of the Wind /SWE (Faraday Cup) data to define 


the orientation of a mean magnetic field, as considered by (Kasper et al. 2008, 2013). Thus 


we can define parallel and perpendicular components of all quantities with respect to the 
direction of the mean field. We will refer to the total proton plasma /3 as the weighed sum 
of its parallel and perpendicular contributions, f3 p = l/3(/? p || + 2/? p j_), where the parallel 
proton plasma /? is defined in terms of the parallel component of the proton thermal velocity 
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and the local Alfven speed, (3 p y = v^ h u/V^. Analogously, the perpendicular part is given 
by /? p j_ = v^JVl For isotropic plasmas the parallel and perpendicular components are 
equivalent: (3 P \\ = /3 P ± = f3 p . We use Wind/ MFI and SWE/Faraday Cup (FC) data to select 
a period with a given proton plasma /?, temperature anisotropy and relative drift speed 
between the protons and the a particles. Once this interval is selected, we parametrize the 
magnetic held fluctuations, which are present there. The SWE/FC data samples are taken 
at 92 seconds, which sets the lowest temporal resolution from all instruments in our analysis. 


The SWE instrument (Ogilvie et al. 1995) provides the proton density, temperature and 


temperature anisotropy, as well as the a particles density and the relative drift speed 
between the two ion species. The electron temperature is obtained from the 3DP instrument 
Lin et al. (1995) on board Wind spacecraft. MFI data ( Lepping et al.||1995 ) provides vector 
magnetic held observations every 0.092 seconds. The different resolution of the instruments 
therefore leads to 10 3 magnetic helds measurements within the 92 seconds period of one 
SWE/FC observation. The analyzed data is from 16 th of January 2008 at 15:19:09 UT and 
represents a long-lasting undisturbed fast wind stream during the recent solar minimum. 
The solar wind speed is 670 km/s. We have identified the time period with (3 P ~ 1 and 
T\\p/T± p ~ 1 and selected the data point closest to (1,1) in the anisotropy-/? phase-space. 
Next, we select a scrolling window of 1000 surrounding points of MFI data and Fourier 
transform them to obtain spectral information about the magnetic huctuations. The data 
analysis has been cross-checked with Morlet wavelet transformations and the same spectral 
properties are obtained. The trace of the power spectral tensor, the magnetic helicity and 
the variance anisotropy were computed using the local magnetic field direction, following 


the procedure explained in Wicks et al. (2012). The normalized relative drift speed for the 


selected interval as obtained from analysis of the data from SWE/Faraday cup particle 


detector is V ap /V\ = 0.44. The proton and the a particle number densities respectively 
are n p = 2.65/cm 3 and n a = 0.09/cm 3 . From the charge neutrality condition the electron 












number density amounts to n e = 2.83/cm 3 . The initial temperature for the a particles has 
not been extracted from the data and for the rest of the modeling we will assume equal 
initial temperature for the two ion species. The amplitude of the mean magnetic field as 
obtained from the 92 seconds average of the Wind/MF\ measurements is |Bo| ~ 4 nT and 
the amplitude of the largest magnetic field fluctuation is about 0.2 nT. Figure [l] shows the 
derived power spectral density of the magnetic held fluctuations in units of nT 2 /Hz as a 
function of the obtained wave frequency given in Hz (in the spacecraft frame) based on 
Fourier decomposition of the MF1 data (Morlet wavelet analysis of the data produces the 
same power spectral slopes.) The level of the instrument noise is over-plotted with a dashed 
line, and the spectra has been fitted by 2 separate power laws /“, where a represents the 
spectral slope. The noise level of the MFI instrument is not always constant nor visible 
in the spectra, but this particular estimate is in-line with the expectations of the Wind 
MFI instrument team (private communications with A. Koval). The original laboratory 
estimation for the instrument suggests a noise-level of 2 x 10 ~ Q nT 2 /Hz. In reality the 
observed instrument noise floor is often slightly lower than that measured in the laboratory. 
For this study we have estimated a value of 8 x 10 ~‘nT 2 /Hz, which is roughly a factor of 3 
lower than what is measured in the lab. The initial lower frequency part of the spectrum, 
associated with large-scale fluctuations shows a / _3 / 2 power law, typical for inertial range 
of a Kraichnan-type MHD turbulence, whereas the higher-frequency (smaller scales) part 
of the spectrum exhibits a steeper power slope a = —3, characteristic for the dissipation 
range. Similar power slopes have been derived for a lower f3 data surrounding the (1,0.33) 
point in anisotropy-/? phase-space (R. Wicks, private communications). In order to study 
direct cascade processes in the present modeling we will initialize the simulations with a 
fraction from the lower-frequency part of the magnetic power spectrum and let the kinetic 
micro-turbulence transfer the wave energy down to the ion scales, where the fluctuations 
can be absorbed due to wave-particle interactions. 


3. Hybrid code — particles initialization 


The details of the 2.5D hybrid code and the numerical schemes used in our model have 


been presented in many previous studies, e.g. Offnan & Vinas (2007); Ofman et al. (2011 


2014); Maneva et al. (2015). Within the present hybrid model the electrons are described 


as a charge-neutralizing massless isothermal fluid, and the ions are treated fully kinetically 
within the particle-in-ccll approach. The code allows for two spatial coordinates, but 
computes all three components of the ion velocities, the current, and the electromagnetic 
fields. The electric and the magnetic fields are computed in Fourier space and periodic 
boundary conditions are assumed for both particles and fields solvers. The normalization 
of the temporal and spatial scales are introduced as follows: the simulation time is given 
in units of inverse proton gyro-frequency 1 and the ion velocities are normalized to the 

Alfven speed, Va = , defined by the electron number density n e , the proton mass 

m p and the magnitude of the homogeneous background magnetic field B 0 . The length of 
the simulation box is given in units of the proton inertial length with equal resolution 
and number of grid points in both special dimensions, so that L x — L y — 384Va/^ p - 
The simulations are performed with 256x256 cells and 127 particles per cell per species; 
thus > 16.6 million particles are used. The time step is a small fraction of the proton 
gyro-period, A t = O.C^flp 1 . The size of each grid cell in x and y direction is the same 
and equals 1.5 proton inertial lengths. The code uses the expanding box model originally 
developed by (Grappin & Velli 1996). Its application to hybrid models was first presented in 


(Liewer et al. 2001) and has lately been widely used by many authors, e.g. Hcllinger et al. 


(2005 

); Hellinger & Travnfcek 

(2006) 

Ofman et al. 

(2011b; 

Moya et al. 

(2012) 

Hcllinger & 


Travnfcek (2013); Maneva et al. (2013); Ofman et al. (2014); Maneva et al. (2015). The 


detailed description of the expanding box model can be found in all of the given references 
and will not be repeated here. For the purpose of the present study we will only introduce 
the solar wind expansion parameter e = Uq/Rq, defined at f 0 — 1/Q p for solar wind outflows 
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with constant velocity U 0 at distance R 0 from the Sun. We load the code with homogeneous 
plasma density and isotropic drifting Maxwellian velocity distribution functions (VDFs) for 
both protons and the minor ions. The initial ion densities correspond to the Wind /SWE 
measurements presented in the previous section, which also provide the initial proton 
temperature, as well as the relative drift speed used to initialize the simulations. The 
electron temperature is obtained from the Wind/ 3DP instrument. Using the electron 
number density as a normalization, the normalized number density for the protons and 
the alpha particles becomes: n p /n e ~ 0.94; n a /n e ~ 0.03. The fluid electron plasma /3 
obtained from the 3DP observations is j3 e = ‘lpQn e k-^T e /B/ = 0.9 The proton plasma (3 is 
defined as the square of the normalized proton thermal speed and is also calculated from 
the Wind /SWE data, see Section |2} We should note that the Alfven speed used here is 
defined in terms of the total electron density n e and the proton mass m p . In this sense 
the definition of the ion plasma j3 for the ion species reads: /3j = m p /rrii2iJ, 0 n e kBTi/B q. 
The plasma /3 for the minor ions is then assumed to be a fraction of the plasma fd for the 
protons, /3 a ~ 0.08; /3 P = 0.33. This selection assures that the two ion species have the same 
temperature at the beginning of the simulations. The selected Wind data interval shows 
equal parallel and perpendicular temperatures for the protons and the same is assumed 
for the a particles. The normalized initial relative drift speed between the two species is 
V ap = 0.44Va- The selected proton plasma-to-cyclotron frequency ratio for the simulations 
is set to ujpi/fl C i = 0.5 x 10 4 , as calculated from the proton density and the magnetic held 
strength provided by the Wind/SWE and MFI instruments. 


4. Construction of the initial turbulent wave spectra 

In the previous section we presented the properties of the hybrid model used for 
this study and the initialization of the particle properties, such as thermal velocities, 
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temperature anisotropies and relative drift speeds. This section describes the theoretical 
reconstruction of an initial Alfven-cyclotron turbulent wave spectra, which resembles the 
Wind /MFI data presented in Figure [lj To construct the 2D oblique wave fronts needed for 
the 2.5D hybrid simulations we first generate a ID broad-band wave spectra and rotate it 
in Fourier space to make it two-dimensional, preserving the energy and the magnetic field 
divergence free. The goal is to start loading the simulations with parallel Alfven-cyclotron 
waves and vary the angle of propagation of the magnetic fluctuations and compare the 
evolution of the related anisotropic ion heating and differential acceleration for the different 
cases. 

Let us assume that the homogeneous magnetic field background is in x direction, B 0 = B 0 cc. 
The general reconstruction of a ID broad-band spectrum of parallel-propagating magnetic 


field fluctuations follows the method developed in Vinas et al. (2014) and Vinas et al. (1984) 

^max 


B (re) = B o + Re 


^ (Ck exp (i<p k ) (y + iz) + exp (i^ k ) (y - iz)) exp (ik\\x) 


L^min 


( 1 ) 


where k min , fc max are the lowest and highest wave-numbers characterizing the wave spectra 
and 
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are the Fourier amplitudes of the left-handed and right-handed modes. The ID spectrum 
describes parallel wave propagation, k = k\\x, with random phases (p k and ^ for each 
mode. In Eq. (|2j) r/ is the maximal spectral wave amplitude, — 1 < a m k < 1 is the reduced 
magnetic hclicity for each mode and a is the spectral index of a prescribed power-law 
power spectral profile. The spectrum conserves the wave energy, polarization, magnetic 


hclicity and the cross-hclicity in configuration space and Fourier space (Vinas et al. 2014) 


This general procedure of generating ID magnetic field fluctuations was applied in a recent 
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2.5D simulation study (Maneva et al. 2015) to construct ID initial broad-band spectra of 


left-hand polarized parallel propagating Alfvcn-cyclotron wave and follow their anisotropic 
nonlinear turbulent cascade. These simulations assumed that the entire simulation box is 
permeated by one-dimensional parallel propagating waves and studied their role for heating 
and acceleration of minor ions in the solar wind. In order to initialize our 2.5D numerical 
simulations with 2D turbulent wave spectra in the present work at y = 0 we start with 
the one-dimensional spectra calculated for all grid points aq along the direction of the 
background magnetic field. Next we transform the resulting ID field B(x) into Fourier 
domain and perform a rotation of angle 9 with respect to the z axis to the obtain 2D 
oblique wave-vectors. The new wave vectors k' and the Fourier component of the complex 
magnetic field B k , acquired through the rotation are given by 


k' = kuX + k' L y = A ;11 cos (9)x + kn sin [9)y . 


( 3 ) 


and 


where 


B' k , = -B k)V sin (9)x + B k , y cos {9)y + B k , z , 


( 4 ) 


( 5 ) 


B k , y = l {C k exp (i<p k ) + Cj? exp (# fc )) 

B k , z = | {C k exp (iip k ) - Cj} exp (iip k )) . 

We apply this procedure to transform the initial ID parallel waves into 2D oblique 
fluctuations at a prescribed angle 9, preserving the energy and the divergence-free property 
of the initial magnetic field V • B ; = 0 (or k 7 • B' = 0 in Fourier domain) and obtain an 
oblique magnetic field Fourier spectrum with power in both parallel and perpendicular 
directions with respect to the background field. In addition, using the new magnetic 
spectrum and wave-vectors we construct the bulk velocities spectra for each ion species, 
namely 

3U s k : y i iUg k^z 


(u/k - U h ) 


( it ( w - ku h ) /n s ) 


B k,y ± iB k,z 


/in 


( 6 ) 









13 


where s is the species index and the +(-) sign represent left (right) hand-polarized waves. 

We repeat the process for each grid row until y = y max , generating N y oblique spectrum 
all with different random phases. To preserve the total energy initially considered we 
must normalize the Fourier coefficients to divide the energy of each perpendicular mode 
to the number of repetitions made in the k' ± direction (by construction the energy is well 
distributed within all the parallel modes k'A. Finally, after the rotation of the entire ID 
Fourier spectrum we obtain oblique two-dimensional spectra of magnetic held fluctuations. 
Simultaneously we rotate and compute the corresponding bulk velocities for each ion 
species as initially given by the parallel dispersion relation, used as our starting case. 
Applying an inverse 2D FFT transform we retrieve the desired oblique 2D waves spectrum 
of Alfven-cyclotron waves with both magnetic and velocity held huctuations, which we load 
in the two-dimensional spatial grid. 


To compute the initial wave frequencies we solve the warm plasma dispersion relation 
for parallel wave propagation in a multi-species homogeneous plasma background. We 
select the left-hand circularly polarized waves and use them to construct the ID velocity 
helds. The determinant of the dispersion tensor in the drifting multi-species magnetized 


plasma is given by the linear Vlasov theory (Davidson & Ogden 1975 Sentman et al. 1981 


Gary et al. 2002) 


det | D | = k 2 c 2 — u 2 — ^ 


CO, 


ps 


co — U s k 
v. s k 


Z(C) + (A S -1)[1 + C^(C)] 


(7) 


The multi-species index s in the summation stands for electrons, protons and He ++ ions. 
The rest of the notations are as follows: oo ps , v s , U s and A s are the plasma frequency, the 
spontaneous velocity, the bulk velocity and the temperature anisotropy for the s-th plasma 
species, A s = Tj_ s /Tjj s . Z(£) is the standard plasma dispersion function with Doppler-shifted 
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frequency-dependend argument 


Z(ff) 


= 5T- 1 / 2 



cu — Ugk T D s 
v s k 


We look for solutions of the above equation for low-frequency waves, with frequencies 
below the lowest ion-cyclotron in the system, in this case the a-cyclotron. For this study 
the electron inertia has been neglected and the electron contribution has been eliminated 
assuming charge neutrality and current conservation. The kinetic plasma dispersion relation 
is used in order to make sure that the selected initial wave spectrum is stable with respect 
to Vlasov theory plasma micro-instabilities. Once the frequencies for the parallel waves 
are obtained, the ID fluctuations of the transverse velocity fields in the initially isotropic 
plasma are calculated from the multi-fluid dispersion relation: 


SV\i = 


oj — Ujk 


SB ± /B 0 


( 8 ) 


1 — (lu — Uik) /D p 

Then the ID velocity fluctuations for each species are rotated to obtain the 2D solution for 
the initial velocity fluctuations in configuration space. The initial wave spectra considered 
here represent two-dimensional energy-preserving Alfvenic fluctuations and more precisely 
forward-propagating left-hand polarized Alfven-cyclotron waves from the lower-frequency 
cc-cyclotron branch of the dispersion relation. The initial wave spectra demonstrated on 
Figure [2] are chosen to resemble the end of the inertial range represented by the spectral 
slope a = —3/2 from the Wind/MFl spectrum, see Figure [lj The generated spectra consist 
of 17 modes and the total amplitude is 5B = 0.2H 0 . For a fixed amplitude and spectral 
slope, we initialize the simulations with three different propagation angles, respectively 
0°,30° and 60°, and study their influence on the minor ions. In the pure parallel case the 
waves have wave numbers k\\ G [0.26, 0.52]D p /Va and frequency range u G [0.22, 0.34]D p . At 
30° the spectrum has k\\ G [0.23, 0.46]D p /Va, oo G [0.21, 0.34]D p and k± G [0.13, 0.26]D p /Va- 
The highly oblique case at 60° has k\\ G [0.13, 0.26]D p /Va, uj G [0.13,0.26]D p and 
k_ |_ G [0.23, 0.46]D p /Va- Figure [2] describes the reconstructed initial magnetic held power 
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spectra as a function of the parallel component of the wave-vector. The left panel shows 
the case of strictly parallel propagating waves, where k = k x x. The right panel shows the 
/^-dependence of the magnetic spectra at fixed perpendicular component k y = 0.13D p /Va 
for the case of oblique wave propagating at 6 — 30°. The same spectral slope a = —3/2 
was used in both cases. Since the initial angle is not too large though, most of the power 
still remains in the parallel direction and the power law retains its slope. This is no longer 
true when highly oblique fluctuations at 9 — 60° are considered. Yet, the magnetic held 
spectra as a function of the full wave-vector remains the same for all angles of propagation 
and recovers the parallel case shown on the left panel of the figure. We should note that 
the initial power spectral law for the magnetic held fluctuation is constructed in /c-space, 
whereas the observational data provides information about the power spectra in frequency 
domain. For the low-frequency MHD part of the Alfvenic spectrum the dispersion relation 
is linear (waves propagate at the Alfven speed) and the two descriptions become identical. 
For the intermediate scales at the beginning of the dispersive part of the spectrum, which 
we address in our study, the spectral slope of the magnetic held huctuations in frequency 
domain which results from k y 7 reconstruction will be hatter with a larger corresponding 
power law exponent a > —1. Figure [3] depicts the 2D power spectral density of the initial 
magnetic held huctuations for the case of oblique Alfven-cyclotron waves propagating at 
6 = 30°. The hgure shows the squared magnitude of the magnetic held huctuations in 
Fourier space as a function of the parallel k x and perpendicular k y components of the 
wave-vector. The initial wave power is restricted in the wave-number space described 
above and presented in Table 1. The following two hgures illustrate the corresponding 
initial magnetic held in configuration space for the case of oblique wave propagation 
at 0 — 30° and 6 = 60°. Figure [4] shows the parallel component of the magnetic held 
huctuations as a function of the parallel x and perpendicular y spatial coordinates, forming 
the simulation box. The huctuations are normalized with respect to the magnitude of the 
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background field B 0 and the spatial coordinates are given in terms of proton inertial length 
di = c/co p i = Va/^ p , where c is the speed of light in vacuum and cu p i is the proton plasma 
frequency. The magnetic field in configuration space for the case of 6 = 30° is given by the 
inverse Fourier transform of the initial fluctuations, presented in Figure [3j The background 
field has not been subtracted, which shifts the magnitude of the fluctuations, forcing them 
to oscillate around the mean normalized value of 1. We should note that in the case of 
strictly parallel wave propagation there are no fluctuation along the background magnetic 
field. The non-zero fluctuations in B x are associated with the presence of oblique waves. 
Figure [5] is a further illustration of the parallel magnetic field oscillations in the case of 
oblique wave propagation at 9 = 30° as an example. The two panel represent a cut of the 
parallel magnetic field component at y — 32di (left panel) and x = 32 di (right panel). As 
these positions the initial waves are still clearly present in both parallel direction along x, 
and perpendicular direction along y regardless the limited range of selected wave-numbers 
for the reconstruction of the initial broad-band spectra, see Figures [3] and [4} 

In the next section, we will present the influence of the initial wave spectra for the 
anisotropic ion heating for the different propagation angles discussed above. We will discuss 
the wave-scattering and the evolution of the ion velocity distribution functions for protons 
and a particles in the electromagnetic field of the waves. We will show the evolution of the 
initial magnetic field and plot the end-stage electromagnetic dispersion relation to analyze 
the scattering of the initial spectra along with the generation of new wave modes. 


5. Results 


In this section we will answer the objectives posed in the introduction, namely to 
compare the relative role of parallel and oblique Alfven-cyclotron waves for the observed 
preferential heating of minor ions in the fast solar wind near the Earth. As a secondary goal 
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of this research, we will investigate the evolution of the stable point /3 P = 0.33; Tj_ p /T\\ p = 1 
in the central domain of the plasma /3 and temperature anisotropy space for the protons, 


where the majority of the ion velocity distribution functions are observed (Matteini et al. 


2007 

Bale et al. 

2009 

Maruca et al. 2011 

2012 Wicks et al. 

2013J Hellinger & Travnfcek 


2014). This point is stable with respect to linear plasma instabilities and it is interesting to 


study its temporal evolution in (j3,T± p /T\\ p ) space in the presence of the ambient solar wind 
turbulence. Figure [6] demonstrates the temporal evolution of the temperature anisotropies 
for the a particles (top panel), protons (middle panel), as well as the relative drift speed 
between the two ion species (bottom panel). The initial plasma state represents the Wind 
measurements explained in Sections 2 and 3. The figure presents a comparison between the 
different simulation cases (see Table 1) with strictly parallel or oblique initial wave spectra 
with propagation angle 6 = 0°, 30° and 60°, respectively. The initial state for all three cases 
is isotropic drifting plasma with V ap = 0.44Va- The slight initial temperature anisotropy 
(mainly in the parallel case) is apparent and is due to the non-thermal contribution in 
the transverse bulk velocity for the ions generated by the initial broad-band wave spectra, 
see Eq. (8). The slight initial deviation from the selected value of the relative drift speed 
in the case of oblique wave propagation is due to the bulk parallel velocity fluctuations 
induced by the oblique waves. The temperature anisotropy of the a particles throughout 
the simulations varies with the propagation angle and is highest for the case of parallel 
wave propagation with Enal (Tj_/Tj|) ~ 1.24, whereas the anisotropy of the protons 

smoothly decreases for all propagation angles considered here. Thus Tj_/Tj| for the minor 
ions increases with 5% for the case of parallel wave propagation, decreases with 5% for 
6 = 30° and is reduced by 20% for the highly oblique waves with 9 = 60°. The anisotropy 
of the protons is reduced for all cases: it decreases with 21% in the parallel case, 11% 
in the case of initial waves at 9 — 30° and 9% for 9 = 60°. These changes in the ion 
temperature anisotropies are related to the preferential heating of the minor ions and 
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the perpendicular cooling for the protons in the decreasing fluctuations of the transverse 
magnetic field. The relative drift speed slowly decreases within the simulation time of 800 
proton gyro-periods, ft decreases with approximately 8% for the case of parallel waves, 
15% for the case of oblique waves propagating at 9 = 30°, and is significantly reduced by 
34% for the case of highly oblique wave propagation at 6 = 60°. We should note that all 
cases above should be stable from the point of view of linear Vlasov instability. Figure [7] 
illustrates the temporal evolution of the parallel (along) and the perpendicular (across 
the ambient magnetic field) components of the ion temperature for protons (top row) and 
a particles (bottom row) throughout the simulation time. The figure is complementary 
to Figure [6] and shows the preferential ion heating for minor ions in the electromagnetic 
field of the non-resonant initial wave spectra. The plots show that the a particles acquire 
high perpendicular temperatures for all cases with an estimated growth rate of 65% to 
72% as we transition from the parallel to the highly oblique initial wave spectra. The 
parallel temperatures for the minor ions are also significantly enhanced and their increment 
throughout the simulations varies from approximately 50% for the initial parallel waves 
through 70% for the slightly oblique case and reaches as high as 120% for the case of highly 
oblique initial wave spectrum at 6 = 60°. The dominant parallel heating in the latter 
case leads to parallel temperature anisotropies for the a particles with Tj_/Tji < 1. The 
intense minor ion heating in both parallel and perpendicular directions can be attributed 
to non-resonant particle scattering and absorption of the initial low-frequency transverse 
fluctuations as well as to Landau and cyclotron damping of resonant a-cyclotron waves 
generated via nonlinear cascade in the course of evolution. The most prominent parallel 
heating of the minor ions in the case of highly oblique waves is due to parallel electric field 
associated with the electromagnetic Alfven-cyclotron waves at oblique wave propagation. 
In the strictly parallel case only longitudinal waves can carry electric field fluctuations and 
the initial spectra of transverse Alfven-cyclotron waves cannot contribute to the electric 
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field and density fluctuations, unless nonlinear wave-wave interactions, such as parametric 
instabilities or more complicated wave couplings take place. Ion heating by parametrically 


unstable Alfven-cyclotron waves has been discussed for example in (Araneda et al. 2009) 


and (Maneva et al. 2014). The protons are less affected by the low-frequency initial 
wave-spectra than the a particles. Their parallel temperature throughout the simulations 
time increases within 5 — 8% and their perpendicular temperature decreases adiabatically 
by 10 — 14% with the decay of the initial transverse magnetic spectra. Figure [8] and Figure [9] 
describe the final stage (at Vt p t = 800) ion velocity distribution functions for protons and 
a particles for the three different angles of propagation of the initial Alfven-cyclotron wave 
spectra. The left panel shows the case of initially strictly parallel waves, the middle panel 
describes oblique propagation at 9 — 30° and the right panel shows the evolved velocity 
distributions for the highly oblique case with 9 = 60°. The bottom rows show the particle 
counts along the parallel direction (solid lines), over-plotted with the best bi-Maxwellian fit 
(dotted lines). The center of all velocity distribution functions is shifted due to the presence 
of relative drift between the a particles and protons. The figures indicate that the final 
distributions for both ion species deviate from the drifting bi-Maxwellian approximation. 
Due to the perpendicular cooling and the parallel heating the velocity distributions for 
the proton species are elongated at all propagation angles with a distinct onset of a beam 
formation in the parallel case. The scattering at the front side (the right hand side) of the 
distribution becomes less prominent in the case of highly oblique wave propagation and 
at 6 — 60° the proton distributions can be sufficiently well-described within the drifting 
bi-Maxwellian approximation. The velocity distribution functions for the a particles 
show strong signatures of particle scattering at the front part of the distributions for all 
propagation angles and at the tail of the distribution functions in the cases of oblique 
wave propagation at 6 = 30° and 6 = 60°. This scattering is caused by the wave-particle 
interactions with the turbulent spectra and relates to the observed perpendicular and 
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parallel heating for the minor ions. The distributions are skewed and throughout the 
simulation time we observe prominent forward beam formation for both protons and the 
a particles in the case of initial parallel wave spectra and a combination of forward and 
inverse a beam formations (referring to the positive/front and negative/back side of the 
velocity distribution function) in the case of oblique wave propagation at 9 = 30° and 
9 = 60°. In the highly oblique case the back scattering by the oblique waves dominates and 
a dominant inverse a beam is formed. Similar beam formations through particle scattering 
by initial parallel wave spectra and monochromatic waves have been reported earlier in 


previous hybrid simulation studies, for example (Maneva et al. 2013), (Maneva et al. 2014) 


and (Maneva et ah 2015). Figure 10 shows the parallel component of the normalized 


magnetic field at the final stage of the simulation, Q p t = 800, for the three different initial 
propagation angles. The scale is given in units of the magnitude of the constant magnetic 
background, whose contribution has not been subtracted from the fluctuations. The top 
panel shows the generated magnetic hied fluctuations in parallel direction for the case 
of initially parallel wave propagation with 9 — 0°. The initial strictly parallel transverse 
wave spectra includes no variations along the direction of the background held and the 
observed small level of huctuations at the hnal stage are entirely due to the generation of 
oblique waves. In the case of oblique wave propagation at 9 — 30° there is a structured 
redistribution of the initial magnetic held huctuations in real space, enhancing the wave 
power in some regions of the box, for instance at x > 300dj and y < lOOd,, see Figure |4j In 
the case of initial wave propagation at 9 = 60° the hnal stage magnetic held also shows a 
structured form different from the shape of the initially constructed huctuations. In both 
cases in the course of evolution the waves start to form more coherent patterns. The total 


amplitude of magnetic held huctuations decreases over time for all three cases. Figure 11 


depicts the evolution of the initial magnetic held power in the Fourier space for the case 
of obliquely propagating initial wave spectra at 6 = 30°. The hgure depicts the strong 
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depletion of the initial wave power with significant direct cascade in perpendicular direction 
particularly visible around k y — k± — 0.36fl p /VA and k y — k± — 0.46h2 p /VA. We should 
note that the initial power is transferred to a wide range of perpendicular wave-numbers, 
which are not initially present in the system (note that the initial range of perpendicular 
wave numbers at 9 — 30° is k± G [0.13, 0.26] fl p /V\). In addition substantial fraction of 
the wave power has cascaded towards strictly perpendicular wave numbers, which implies 
the generation of dispersionless wave modes propagating at 9 = 90°. As visible in the next 
figures this process persists for the case of initial spectra at 9 = 60° where the wave energy 
is transferred to even smaller transverse scales with higher perpendicular wave-numbers. 
Apart from the perpendicular cascade the spectra also shows both direct and inverse 
cascade in parallel direction. The direct cascade happens at earlier stage before the end of 
the simulations and generates resonant modes, which are quickly absorbed by the minor 
ions. Figures p~2| through p~4| provide information about the nonlinear evolution of the initial 
wave-spectra and the dispersion of the turbulence generated waves in the course of the 
simulations. The simulation dispersion relation is computed from the power of the Fourier 
transformed magnetic field fluctuations in space and time. To compute the dispersion 
in parallel direction, the Fourier transformation is taken at a fixed position, from the 
vertical domain of the initial wave spectra. The same procedure is used in the computation 
of the wave dispersion in perpendicular direction. To compute the dispersion we have 
used the temporal information throughout the end of the simulations. The magnetic 


held power along the vertical line perpendicular to in the bottom panel of Figure 12 


is due to the strong parallel waves, plotted in the top panel. The horizontal lines with 
to < 0.5h2 p may relate to ion Bernstein modes for the a particles, while higher-frequency 
proton-Bernstein modes are highly damped. The figure indicates the generation of strong 
forward propagating fast modes, as well as strong backward propagating Alfven-cyclotron 
modes at wave-numbers range spreading beyond the wave-numbers of the initial parallel 
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wave spectrum. Low-frequency backward propagating Alfven waves are undamped solutions 
of the parallel linear Vlasov dispersion relation when sub-Alfvenic drifts are present. Here 
they are excited through nonlinear wave-wave interactions in the presence of the initial 
wave spectra and the induced ion temperature anisotropies. The generation of these 
background propagating waves continues at oblique wave propagation, but they are limited 
within lower wave-number range and at 6 = 60° their signature is overcome by strong 
backward-propagating fast modes, as visible from the next figures. The forward propagating 
fast modes are less prominent for the case of 9 = 30° and are overcome by backward 
propagating fast modes when 6 = 60°. The zero-frequency lines at the bottom panels in 
Figure [13] and Figure [14] describe non-propagating entropy waves, generated by the density 
and pressure fluctuations. In an isotropic and slightly anisotropic plasma with sub-Alfvenic 
drifts within the specific wave-number range linear instability theory predicts that these 
waves should be stable or fully damped. The presence of the relative drift speed reduces the 
damping of the slow modes and the entropy waves and the pressure fluctuations associated 
with the initial wave spectra further enhance the entropy waves wave-power. Thus the 
positive growth rate of the entropy waves is a nonlinear effect, which can be attributed to 
the presence of the initial Alfven-cyclotron wave spectra. Looking at the complementary 
spectral plot shown on Figure [IT] and the explanations thereafter, the zero-frequency modes 
can be attributed to the generation of strictly perpendicular Alfven waves. 


6. Discussion and Conclusions 

Parallel and obliquely propagating ICWs are ubiquitous in the solar wind plasma 
and have been observed at various heliocentric distances. Furthermore, differentially 
streaming non-thermal ions with beams in their velocity distribution functions and different 
temperatures along and across the direction of the ambient magnetic held are common 
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features of the fast solar wind streams. Yet, the correlations between the existing plasma 
waves and the particle properties in the multi-species solar wind are still poorly understood, 
and further theoretical and simulation studies directly related to observations are needed. 
In this paper we present the results from observations-driven 2.5D hybrid simulation study, 
as a first attempt to distinguish the relative role of parallel versus oblique wave propagation 
for the generation and evolution of ion temperature anisotropies and relative drifts in the 
fast solar wind streams. The initial state is based on Wind data from 16 th of January 
2008, 15:19:09 UT time in the undisturbed wind at 1AU, and we study the evolution of 
the expanding solar wind plasma parcel, which cannot be followed by the spacecraft. We 
have investigated the evolution of a stable region within the velocity distribution functions 
of proton and a particles in anisotropy-/? space. The initial isotropic distributions of the 
minor ion species interact with the wave spectra and acquire slightly higher perpendicular 
temperatures. In contrast the protons cool down, as there is no wave energy at the 
proton scales. Our findings show that the minor ions are preferentially heated by the 
magnetic fluctuations for all initial wave spectra considered here, where protons are prone 
to perpendicular cooling, as expected by the double-adiabatic expansion model predictions. 
At the final stage of the simulations the strength of the magnetic fluctuations decreases 
as the magnetic power is converted into thermal energy for the minor ions. The main 
assumption in this work is the type of initial wave spectra considered here. Since we are 
interested in the properties of the minor ions and try to couple the physics at fluid and 
ion scales, for the initial broad-band wave spectra we have chosen low-frequency left-hand 
circularly polarized Alfven-cyclotron waves. We exclude the effect of the high-frequency 
(much higher than the proton-cyclotron frequency) whistler and KAWs waves, which 
interact mainly with high-speed electrons and proton beams. Instead our study is focused 
on the direct cascade of wave power from the intermediate fluid regime down to the kinetic 
ion scales and neglect the effects of inverse cascade at scales much smaller than the proton 



24 


skin depth or their possible contributions to the energy transfer to the ion scales. Yet, apart 
from the direct cascade towards the smaller ion scales, our simulations do show signatures 
of inverse cascade towards larger fluid scales with generation of strictly perpendicular 
waves as clear from the magnetic held power spectra plot at the end of the simulations. 

In this study we have performed 2.5D hybrid simulations to investigate the minor ion 
heating by broad-band spectra of oblique Alfven-cyclotron waves and the evolution of an 
initial relative drift speed between the protons and the minor ions. We have compared the 
ion heating corresponding to different angles of propagation (0°,30° and 60°) and found 
prominent perpendicular heating for the minor ions at all propagation angles. For the 
selected spectral range of the initial wave spectra for the simulations we find that the a 
particles are most prominently heated in both parallel and perpendicular direction by the 
highly oblique Alfven-cyclotron waves propagating at 6 = 60°. The parallel heating for 
both species increases with the propagation angle and the parallel temperature of the minor 
ions in the highly oblique case prevails over the perpendicular one, so that they acquire 
parallel temperature anisotropy with Tj_/Tj| < 1. Since the perpendicular heating of the 
minor ions remains at similar level for all propagation angles the minor ions have highest 
anisotropy for the case of initially parallel propagating waves. For all cases considered here 
the initially isotropic protons experience moderate cooling in perpendicular direction, and 
slight heating in parallel direction, the net effect of which leads to the generation of parallel 
temperature anisotropies. Contrary to the linear instability theory predictions the initial 
relative drift speed V ap = 0.44Va decreases in time throughout the simulations. It remains 
quasi-conserved in the parallel case and gets further reduced in the oblique cases with a 
magnitude change of over 30% when 6 = 60°. As part of this study we have also looked at 
the effect of a gradual solar wind expansion with et 0 = 10~ 5 . As expected by construction, 
the expansion leads to perpendicular cooling for both species, however considering the 
realistically small value of the expansion parameter it does not have significant effect on the 



25 


ion heating within the simulation times considered here. The current study is concentrated 
on warm plasmas with f3 « 0.33, but this result is expected to hold even stronger when 
applied to the lower plasma f3 conditions in coronal holes and the fast solar wind near the 
Sun, where /3 ~ 10~ 2 . Next, we investigate the power spectra of the magnetic fluctuations. 
We find limited turbulent cascade in the perpendicular direction for the parallel case, 
accompanied by a strong inverse cascade in parallel direction towards smaller wave-numbers. 
For the oblique cases we find prominent additional oblique mode generation in a wider range 
of wave-numbers than the initial values, as well as the generation of strictly perpendicular 
waves. Right-hand polarized fast waves are being generated in all simulation cases with 
strong forward propagating modes in the parallel case and strong backward-propagating 
modes in the highly oblique case. Strong forward-propagating ion-cyclotron modes are 
excited in the parallel case, as they are expected normal modes in drifting two ion-species 
systems. In the cases of initial oblique wave propagation the initial Alfven-cyclotron wave 
spectrum is strongly damped and generation of entropy modes/perpendicularly propagating 
Alfven waves is observed. We should note that the total magnetic energy is gradually 
depleted throughout the simulation time for all three initial angles of wave propagation 
and it decreases approximately by 70% from the initial value. Simultaneously the minor 
ions are continuously heated on average by 69% in perpendicular direction. This analysis 
suggests that major part of the minor ion heating is due to the external wave energy sources 
and can not be simply attributed to redistribution of kinetic energy between the protons 
and the a particles. Further source for the ion heating can be a nonlinear magnetosonic 
instability induced by the initial drift-speed in the presence of the wave spectra. The value 
of the initial relative drift is within the stable regime for linear plasma instabilities, yet it 
decreases up to 30% in the highly oblique case. This can contribute to the energy source 
of the initial wave spectra and provide an additional energy source needed for the strong 
parallel minor ion heating with more than 120% increase of the parallel temperature for the 



26 


minor ions. Note that the relative decrease of the proton temperature anisotropy, which 
varies between 13% and 20%, is due to a combination of parallel heating and perpendicular 
cooling. The drop in proton anisotropy is most prominent for the parallel case and weakest 
for the highly-oblique case, whereas the minor ion heating is strongest for the highly oblique 
case. This is a clear indication that the observed minor ion heating is related to the initial 
wave-spectra, rather than energy exchange between the two ion species. 

Looking into the velocity distribution functions we find strong signatures of wave scattering 
with prominent forward or inverse minor ion beams being formed by Vt p t = 800. The 
proton scattering by the oblique waves generates leads to no clear beam formations. This 
can support the conclusion that stronger proton beams are formed by trapping in the 
ion-acoustic waves potential generated in the nonlinear evolution of the initial wave-spectra 
or relate to the fact that it is more difficult to accelerate protons by low-frequency waves 
in the finite-/? plasma considered here in comparison to the lower-/? regimes close to the 
Sun, where the wave amplitudes are also expected to be stronger. We should note that 
for the limited wave-number range considered in our study the initial spectral slope of the 
fluctuations does not play an important role for the nonlinear wave energy transfer and the 
related ion heating/cooling and acceleration or deceleration. Simulation tests with more 
gradual initial magnetic field power slope, for instance a = —1, show similar results to the 
ones presented here. This work does not include post-simulation wave polarization analysis, 
however based on previous experiments and the presented wave dispersion analysis in 
Fourier space we believe that right-hand polarized waves are being generated in the course 
of evolution and the final state consists of a mixture of left- and right-handed waves. In a 
next improvement of the present analysis we will broaden the spectral range of the initial 
fluctuations and include different polarization states. As a future perspective we plan to 
implement a two-dimensional spectrum of a general type, avoiding the rotation of parallel 
spectra. We would also consider low-frequency right-hand polarized kinetic Alfven waves 
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and whistlers, and investigate the amount of energy they can deposit to the ions before 
ultimately dissipating at the electron scales. 
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Table 1: Initial parameters for the different numerical simulation runs. The table presents the simulation 
cases as a function of the angle of propagation of the initial Alfven-cyclotron wave spectra, 6 , the initial 
relative drift speed, V ap , the range of initial wave frequencies u> o, parallel, fey, and perpendicular, k±, wave- 
numbers. _ 


Case # 

e 

V ap [V A ] 

“0 [Op] 

fc|l [Op/V A ] 

[Op/V A ] 

1 

0° 

0.44 

[0.23-0.34] 

[0.26, 0.52] 

0 

2 

30° 

0.44 

[0.21-0.32] 

[0.23, 0.46] 

[0.13, 0.26] 

3 

60° 

0.44 

[0.13-0.26] 

[0.13, 0.26] 

[0.23, 0.46] 
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Fig. 1.— Magnetic field power spectrum based on Wind/MFl data fitted with two different 
power laws. The data is from 16th of January 2008, 15:19:09 UT. 



Fig. 2. — Initial magnetic field power spectra as a function of the parallel wave-number for 
the case of parallel (left panel) and oblique (right panel) wave propagation at 6 = 30°. The 
oblique power spectra is calculated at k y = 0.13O p /Va, but the spectral slope is the same 
for the entire interval of initial perpendicular wave-numbers k y € [0.13,0.26]f2 p /VA- 
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Fig. 3.— Power spectral density for the initial two-dimensional magnetic field fluctuations 
as a function of the parallel, k x , and perpendicular, k y , wave-numbers for the case of oblique 
wave propagation at 9 = 30°. The magnetic field power is concentrated within the selected 
wave-numbers area presented in Table 1. 










30 



Fig. 4.— Parallel component of the two-dimensional magnetic field in normalized units 
in real space. The field corresponds to oblique propagating waves at 6 = 30° (left) and 
6 = 60° (right). The spatial coordinated are given in units of proton inertial length. The 
fluctuations are superimposed on the constant magnetic field background, whose amplitude is 
used for normalization. The fluctuations introduce inhomogeneity in the uniform background 
magnetic field. 




Fig. 5.— Parallel component of the two-dimensional magnetic held in normalized units 
as a function of the parallel (left) and perpendicular (right) spatial coordinates, x and y. 
The fluctuations correspond to obliquely propagating waves at 6 — 30°. The plots are 
complementary to the left panel of Fig. [4j 
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Fig. 6.— Temporal evolution of the ion temperature anisotropies Tjj/Tju and relative drift 
speed for the case of parallel (solid red line) and oblique wave propagation at 6 — 30° (blue 
dashed line) and 0 = 60° (brown dash-dotted line). 
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Fig. 7.— Evolution of the parallel (dash-dotted blue lines) and perpendicular (solid red) ion 
temperature components for a particles (top row) and protons (lower row). The different 
panels show the different perpendicular heating for the minor ions induced by the 3 different 
initial wave spectra of parallel (first column) and oblique waves at 6 = 30° (second column) 
and at 6 = 60° (third column). 
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Fig. 8.— Snapshots with the velocity distribution functions for protons for initially parallel 
waves (left panel), and oblique waves at 6 = 30° (middle panel) and at 6 = 60° (right panel). 
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Fig. 9.— Snapshots with the velocity distribution functions for a particles for initially 


parallel (left panel) and oblique waves at 6 = 30° (middle panel) and at 6 = 60° (right 
panel). 
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Fig. 10.— 2D magnetic field fluctuations along the background field at Vt v t = 800 as a 
function of normalized coordinates in configuration space. The top panel relates to oblique 
fluctuation formed from the initially parallel wave spectra. The middle panel shows the 
end-state magnetic field variations for the case of oblique waves propagating at 9 = 30 deg. 
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Fig. 11.— Power spectra of the evolved magnetic field fluctuations at the final stage of the 
simulations at £} p t = 800. The figure shows the square of the magnetic field fluctuations as 
a function of normalized parallel and perpendicular wave-numbers for the case of initially 
obliquely propagating waves at 6 = 30°. 
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Fig. 12.— Power spectra of the magnetic field fluctuations in Fourier space at the final 
stage of the simulations at Q p t = 800. The figure shows the square of the magnetic held 
fluctuations in logarithmic scale as a function of normalized frequency, parallel (top panel) 
and perpendicular (bottom panel) wave-numbers for the case of initially parallel wave prop¬ 
agation. 
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Fig. 13.— Power spectra of the magnetic field fluctuations similar to Fig. 12 but for the 


case of oblique waves propagating at 9 = 30' 
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The same spectral plot for the magnetic field fluctuations as on Fig. 12 and 


Fig. 

Fig. 


but for the case of initial oblique wave propagation at 9 = 60°. 
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